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SUMMARY 


A brief study of the Zimmerman flutter -margin criterion has been made by apply- 
ing it to data obtained from a wind-tunnel model. The sensitivity of the flutter- 
margin parameter was explored with a parametric trend study and by calculation of the 
derivatives with respect to the input frequency and damping parameters. The crite- 
rion is simple in concept and application, and it serves as a good flutter-onset 
predictor because it gives a nearly linear variation with dynamic pressure. However, 
accurate values of both frequency and damping of both modes involved in flutter are 
required for reliable flutter-onset prediction. The simplified version using only 
frequencies gave a highly nonconservative flutter onset in one case and should not be 
used in general. Normalizing the flutter-margin parameter by the wind-off values of 
the simplified flutter-margin parameter yields a parameter of the order of 1 .0 to 0, 
and this is the recommended procedure. 


INTRODUCTION 

One difficult aspect of flutter testing is that a flutter mode can suddenly 
become unstable with only a small increase in dynamic pressure. Furthermore, there 
may be little or no indication of the approach of the instability in a plot of modal 
damping against dynamic pressure. Although such a sudden onset of instability is . not 
always the case, it occurs in practice frequently enough that flutter prediction 
based on the projection of modal damping may not be adequate. The development of 
techniques for determining the nearness of flutter onset has received considerable 
attention. (See refs. 1 to 8, for example.) These flutter-prediction techniques, 
which vary considerably in form and complexity, extend from projections of inverse 
peak amplitude of modal response to ambient excitation (ref. 7, for example) to pro- 
jections based on the differential equations of motion using coefficients determined 
from the dynamic measurements of forced response (ref. 3) . The technique to be used 
in practice depends on such factors as the availability of input forcing, instrumen- 
tation, the amount of computation and analysis required, and whether an on-line or 
off-line technique is required. In many cases, no single technique is sufficient and 
the use of several techniques is required for reliable flutter-onset prediction. The 
purpose of this paper is to present the results of a brief study of one technique 
that was previously developed (ref. 2) and that is commonly referred to as the 
Zimmerman flutter criterion. 

The Zimmerman flutter criterion is based on considering the value of Routh's 
discriminant for the characteristic equation of a two-degree-of-f reedom system as a 
flutter-margin parameter. It can be expressed in terms of the measured frequency and 
damping of the two modes which are input as parameters. Although this technique was 
developed over 15 years ago, it has not been widely used because it requires measured 
values of the frequencies and dampings of both modes involved in the flutter con- 
dition. Typically, the damping of one of the modes increases so rapidly as flutter 
is approached that it becomes difficult to measure and is, thus, unavailable. With 
modern testing and data-analysis techniques, this difficulty should be alleviated and 
the criterion should be of more interest than in the past. 

The flutter-margin criterion and its analytical development are described in 
this report. The criterion is applied to a wind-tunnel flutter model for which some 



analytical and experimental data are available. (See ref. 9 .) The sensitivity of 
the flutter margin to the modal frequencies and dampings is evaluated with a para- 
metric trend study and with the development of derivatives of the f lutter-marqin 
parameter. 


SYMBOLS 

Aj coefficients of characteristic equations, where j = 0,1, 2 , 3 

c chord 

F flutter-margin parameter normalized by simplified value of F s ^ 

F flutter-margin parameter 

F simplified flutter-margin parameter normalized by p _ 

S S f u 

F simplified flutter-margin parameter 

s 

F s 0 simplified f lutter-marqin parameter evaluated at q = 0 
f frequency, o)/ 2 -rc, Hz 

i 

M free-stream Mach number 

q free-stream dynamic pressure, psf 

P real part of complex characteristic root (decay rate or damping), 1/sec 

X complex characteristic root, |3 + iw 

w imaginary part of complex characteristic root (frequency), rad/sec 

Subscript: 

i modal index, i = 1,2 


ANALYSIS 

Zimmerman Flutter-Margin Criterion 

The stability of a two-degree-of-f reedom flutter system with quasi-static aero- 
dynamics is given by the quartic characteristic equation 

X 4 + A3X 3 + A 2 X 2 + AX + A q = 0 ( 1 ) 


2 



where the values of 


A are real 


The roots X are usually two complex pairs 


N,2 ' 51 l “l 


X = P + iut 
3,4 k 2 2 


The boundary for neutral stability is given by Routh's criterion 


A 


2 





= 0 


( 2 ) 


For a stable condition, this equation is not satisfied and the left-hand side of 
equation (2) yields a positive number F, which is defined as the Zimmerman flutter- 
margin parameter and is generally called the Zimmerman flutter criterion. The values 
of A can be expressed in terms of p and w by expanding the factored equation 

( X - P 1 + io^MX - p 1 - iw 1 )(\ - p 2 + ito 2 )(X - P 2 - iw 2 ) = 0 (3) 


and comparing the terms with those of equation (1). Then, F can be expressed in 
terms of p and w as follows: 



Note that 


1 




0 


if either P 1 or P 2 = 0 and, consequently, F = 0. Equation (4) is one of the 
many algebraic forms of the Zimmerman criterion (eg. (A21) of ref. 2) that can be 
obtained by manipulation of equation (A20) of reference 2. 
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A simplified version of equation (4) can be derived by using only static aero- 
dynamics, where all modes are undamped up to the speed for frequency coalescence. 
Letting (3^ = |3 0 in equation (4) gives the result 



where the subscript s denotes the static aerodynamic value. 

The values calculated for F in equation (4) tend to be large and of the order 
of w 2 4 . A convenient way of normalizing f is by substituting the value of F 
at q = 0. Thus, by letting 



and 

F = F / F ~s,0 
the result is 



which gives a parameter that is normally in the range from 1 ."0 to 0. Normalizing 

F in the same manner qives 
s 



( 6 ) 


Note that F s is always positive and is zero only when = co-j • 

For flutter-onset prediction, F is calculated from measured values of 
P and to for several discrete values of dynamic pressure at constant Mach number. 
A least-squares fit of a straight line can then be used to project to flutter onset 
(where F = 0). The fit and the onset prediction can be updated as further points 
are measured. 
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It can be shown that for quasi-static aerodynamics and without structural damp- 
ing, F is a quadratic equation with dynamic pressure q. (See ref. 2.) Use can 
be made of the quadratic relation in order to project to flutter. However, with 
unsteady aerodynamics the relation is more complex than quadratic and, thus, this 
relation is not used herein. If aerodynamic lags of higher order than quasi-static 
are considered (such as exponential lags for the indicial response), the quartic 
equation (1) becomes a polynomial of higher degree. 

Derivatives of F 

The sensitivity of F to all the input values of ^ can be determined by 
differentiating equation (5) and calculating dF/df^ . The resulting equations are 


dF 

da) 


4 M 2 “l 



(7a) 


dF 

dco 2 


4 Pl^2 aJ 2 1 



(7b) 


dF 

dp. 


M 2 * 5 ,'?, + p 2 > 



(7c) 



P, [2P 2 (P, 



(7d) 


where 


2^2 


„ >2 _/■ 2 2> ^2 °1 / 2 2 ^ 
W = ^1 + ^2 + + w 2 ' + 3 l t °2 " ) 

(Pi + p 2 ) 


( 8 ) 
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and dF/d^ are evaluated 


The derivatives dF/dto-j 
and f 2 , in hertz, by 


respectively, in terms of f 1 


dF 

df 


1 


2n 


dF 

d 


(9a) 




(9b) 


Similar results for F g are as follows: 


dF_ 


du. 



(10a) 


dF 


s 


dw 2 



(10b) 


2 

Alternative derivatives with respect to w may be more appropriate since F s 
depends only on the squares of and w 2 and are, respectively. 


dF 


s 




(11a) 





(11b) 


RESULTS AND DISCUSSION 
Description of Model and Tests 

The sample case used for application of the flutter-margin criterion is the 
wind-tunnel model of references 9 and 10. The planform is shown in figure 1. The 
wing was of conventional flutter-model construction with a spar and pods and was 
cantilevered from a sidewall mount. It was flutter tested in the Lanqley Transonic 
Dynamics Tunnel by varying the test-medium density (Freon 1 ) at a constant Mach number 
for Mach numbers of 0.6 and 0.9. Experimental flutter points are given in refer- 


1 


Freon: 


Registered trademark of E. I. du Pont de Nemours & Co., Inc. 
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ences 9 and 10. Calculations for subcritical values of density at M = 0.6 and 0.9 
are given in reference 9 and are used herein. The subcritical calculations are based 
on fitting doublet lattice aerodynamic forces for harmonic motion with a ratio of 
polynomials in reduced frequency, and then on using the fitted function to generalize 
for transient motions. (See refs. 9 and 10.) In addition, some limited, previously 
unpublished, experimental data at subcritical densities were available. The frequen- 
cies of the first two modes were determined from on-line spectral measurements of 
accelerometer response to wind-tunnel turbulence. Ensemble-averaged spectra in the 
range from 0 to 50 Hz were taken by using 20 ensemble averages (102 sec of data). 

The frequencies of the first two modes were determined from the peaks. Damping and 
frequency values for the first mode were also determined from fast frequency sweeps 
of the control surface (from 2 to 20 Hz in 20 sec) for subcritical conditions at 
M = 0.9. The records were analyzed by using the methods and equipment described in 
reference 11. No information on the damping of the second mode was available. The 
flutter data (refs. 9 and 10) indicate that the bottom of the transonic dip is above 
M = 0.9. The first two mode frequencies calculated by a finite-element method 
(ref. 9) are 5.23 Hz and 19.13 Hz. The corresponding measured values are 5.2 Hz 
and 19.2 Hz. The first mode is first wing bending, and the second mode is a com- 
bined bending and torsion mode. Although references 9 and 10 describe active 
flutter-suppression tests, only the data for the wing with the flutter-suppression 
system off are considered herein. 


Flutter-Margin Calculations 

The calculated frequencies and decay rates for the first two modes (from ref. 9) 
and the corresponding values of F and F g are shown in figures 2 and 3 for Mach 
numbers of 0.6 and 0.9, respectively. The corresponding experimental flutter points 
(figs. 2 and 3) are in good agreement with the calculated values. For both Mach 
numbers, the variation of F with q produces nearly a straight line, except near 
q = 0, and, thus, would serve as a much better projector to flutter than modal decay 
rate or damping. For example, for M = 0.6, damping of the flutter mode changes 
curvature and slope and then becomes unstable between q = 130 and 155 psf, whereas 
F is nearly a straight line for q > 80 psf. It might also be noted that the calcu- 
lated frequency and damping values are based on solutions using 10 calculated vibra- 
tion modes and unsteady aerodynamics. The influence of the vibration modes above the 
second mode on the calculated values of F is apparently small in this case. 

For M = 0.6 (fig. 2(c)), F g is quite close to F, except near q = 0. For 
this case, the modal frequencies are, in fact, near coalescence at flutter. However, 
F g does not go to zero but remains a small positive number since the modal frequen- 
cies are not equal at flutter as assumed for F g . For M = 0.9 (fig. 3(c)), F g 
is not close to F because the modal frequencies are not near coalescence, even at 
flutter. The use of F g as a predictor would overestimate the flutter dynamic pres- 
sure by 30 percent. This unconservative estimate indicates that, in general, F g 
should not be used for flutter-onset prediction. 


Derivatives of Flutter-Margin Parameter 

The derivatives of F with respect to the input parameters , |3 2 r f ^ , 

and f 2 are presented in figures 4 and 5. These derivatives are calculated from 
equations (7) to (9) for the two Mach numbers and the dynamic pressures of the sample 
case by using the frequency and damping (decay-rate) values from figures 2 and 3. 
These derivatives indicate the local sensitivity of F to small changes in the input 
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parameters (such as measurement error) and indicate the emphasis that should be given 
to each measurement. An extreme sensitivity would indicate a difficulty with the 
method. 

At the flutter point, all the derivatives, except dF/dp^ , are zero because p^ 
is zero. Thus, p i is the only significant parameter near flutter for this case 
where flutter occurs in the first mode. As previously indicated in the "Analysis" 
section, F = 0 if either p 1 or p 2 = 0; thus, the low damping of the flutter 
mode determines F at or near the flutter point. The primary concern, however, is 
not at the flutter point but with the sensitivity of F at subcritical speeds for 
use in projecting to flutter. The derivatives of dF/dp 1 and dF/dp 2 are very 
large near q = 0 where both p 1 and p 2 are small (essentially structural 
damping) . This result indicates that the large differences in F and F s near 
q = 0 (figs. 2 and 3) result from the inclusion of damping in F. For this case, 
the derivatives dF/dp^ and dF/dp 2 are large where the frequencies of the two 
modes are well separated, and they are small when the frequencies are near coales- 
cence. If a measurement error of 1.0 rad/sec is assumed for p^, the change in F 
for 60 psf < q < 100 psf would be relatively small for M = 0.6 (fig. 4), but it 
would be of the order of 0.1 to 0.2 for M = 0.9. Although dF/dp 2 is smaller than 
dF/dp.|, the measurement error would be expected to be correspondingly larger where 
P 2 is large and its measurement more difficult. Thus, F would also be signif- 
icantly influenced by p 2 away from the flutter point. 

The derivative dF/df 2 is generally much larger than the derivative dF/df 1 
(figs. 4 and 5) and differs in sign; this indicates that F is more sensitive to 
f 2 than f -| . The derivatives of F g , the simplified criterion which is given by 
equation (10), indicate that dF s /df 2 should be larger than dF s /df 1 by the ratio 
of the two frequencies and it should be of opposite sign. Similar results for F 
are apparent in figures 4 and 5. The value of dF/df., is nearly a constant at small 
values of g and is small when considering a reasonable measurement error for f 1 
of 0.25 Hz. The value of dF/df 2 is larger at small values of q where the fre- 
quencies are well separated, and it decreases to become of nearly the same order of 
magnitude as dF/df^ as the two frequencies approach the flutter condition. 

Overall, the sensitivities indicated by figures 4 and 5 are reasonable; and any 
scatter in F would be expected to be of the order of magnitude of any scatter in 
p and to, except for the damping parameters near q = 0. 


Parametric Trend Study 

The expression for F (eq. (5)) is a simple algebraic function of the four 
parameters p 1 , p 2 , io.j , and to 2 . The influence of each parameter can be evaluated 

by varying each parameter with the other three held constant. The results, which are 
presented in figure 6, complement the local sensitivities given by the derivatives 
discussed in the previous section by indicating the change in F for finite changes 
in the parameters. 

The variation of F with p^ for several values of p 2 is shown in fig- 
ure 6(a) for values of frequency at q = 0. The frequencies for this case are well 
separated and, as previously indicated, F is very sensitive to p 1 for small 
initial values of p^ or p 2 * Moderate sensitivity is still indicated (by the 
slope of the curves) for moderate values of both p 1 and p 2 , and less sensitivity 
is indicated for large values of either p^ or p 2> Similar trends are noted for 
varying p 2 for constant values of p,j. (See fig. 6(b).) 
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The variation of F with f 1 for constant values of f 2 is shown in fig- 
ure 6(c), and the variation of F with f 2 for constant values of f 1 is shown in 
figure 6(d). For these cases, F is normalized by F g Q with f 1 = 5.23 Hz and 
f 2 = 19.13 Hz. The value of F, therefore, varies considerably as f.| and f 2 are 
varied. For this case with small values of p.j and p 2 , F is sensitive to f 1 
only when f^ approaches f 2 (fig. 6(c)), where F becomes nearly zero. Similar 
trends are noted as f 2 is varied (fig. 6(d)), but there is considerable sensitivity 
to f 2 (large slope) for larger values of f 2 as previously indicated by the 
derivatives . 


Calculations Using Limited Experimental Data 

A limited amount of subcritical frequency and damping data for the first two 
modes were available from the experiment. In this section, flutter-onset prediction 
is examined by using a combination of the available experimental data and the analyt- 
ical results. 

The frequencies of the first two modes from ensemble-averaged spectral measure- 
ments are compared in figure 7 with the calculated values from figures 2 and 3. 

There is good overall agreement, although there are some differences in the values of 
the first-mode frequency, and there is some scatter in the measured frequencies for 
the second mode. The corresponding values of F calculated by using the measured 
frequencies and the calculated damping values (from figs. 2 and 3) are shown in fig- 
ure 8. The use of the measured frequencies results in little chanqe in F except 
for the lower dynamic pressures. As previously indicated, F t 0 as p 1 -> 0 and p 
determines the value of F near the flutter point. The measured frequencies, thus, 
have little influence near the flutter point. At the lower values of dynamic 
pressure, using the measured frequencies for this case would lead to an underpre- 
diction of the flutter point. Thus, using the measured frequencies does not signifi- 
cantly improve flutter-onset prediction. 

The values of frequency and damping of the first mode as determined from fast 
frequency sweeps are compared in figure 9 with the calculated values from figure 3 
for M = 0.9. The measured frequencies are in qood agreement with those from the 
ensemble-averaged spectral measurements. (Compare figs. 7 and 9.) The measured 
damping values show the same trend as the calculations and are in good agreement with 
the calculated values at low dynamic pressures. They differ near the flutter point 
because the calculations predict a somewhat lower dynamic pressure at flutter. The 
corresponding values of F calculated by using the measured results for the first 
mode and the calculated results for the second mode are presented in figure 10. They 
are also compared with the results obtained from using calculated values of both the 
first and second modes. The experimental flutter point is closely defined by the 
values of F (fig. 10), whereas the measured values of do not define the 

experimental flutter point nearly as well (fig. 9) . Flutter-onset prediction from 
the values of F for q < 80 psf would be impossible, however, as a result of the 
scatter in the measured values. (See fig. 10.) 

It appears that reliable onset prediction requires accurate measurement of the 
frequency and damping of both modes. Although mixing some calculated and measured 
results was helpful in one case, the general use of this technique is questionable. 



General Comments on the Technique 


One of the factors discussed in the previous sections is the sensitivity of the 
results to the input damping and frequencies. Overall, the sensitivity was reason- 
able and the scatter in F would be of the same order as that in the individual 
damping and frequency values used to compute F. However, the plot of F against 
q is nearly linear and would, thus, be less sensitive to scatter than would the 
complex curves that are obtained by plotting damping against q. Some improvement 
may also be gained by considering projections based on a quadratic function for F 
plotted against q as can be seen in figure 2(c), where F plotted against q is 
somewhat nonlinear. It is preferable to use linear extrapolations wherever possible 
and to limit nonlinear projections to a very small range of projection. Thus, 
caution should be exercised if the quadratic projection is used. 

The brief results of this report have also indicated that further emphasis should 
be given to the accurate measurement of the frequency and damping of the second mode 
in addition to the flutter mode. In general, this is a difficult task and may require 
special excitation techniques when the damping of the second mode is large. 

The application cited here and the development of the criterion presented herein 
have considered the case of flutter involving only two aeroelastic modes. It might 
be noted, however, that at least one case has been reported in the literature 
(ref. 4) where the projection near flutter was quite nonlinear, possibly as a result 
of the involvement of more than two modes in the flutter condition. The equivalent 
Routh's discriminant can be readily derived for more than two modes, but its appli- 
cation to flutter data has not yet been reported. In addition, other methods have 
recently been developed (ref. 8) that may complement the use of the Zimmerman 
criterion for multimodal situations. 

The result of the application of the Zimmerman criterion to wind-tunnel and 
flight flutter testing is apparent. Other applications may be for the projection to 
flutter from the output of time-dependent, transonic, aeroelastic analysis codes, 
such as that of reference 12. 


CONCLUSIONS 

A brief study of the Zimmerman flutter-margin criterion has been made by 
applying it to data obtained from a wind-tunnel model. The sensitivity of the 
flutter-margin parameter was explored with a parametric trend study and by 
calculation of the derivatives with respect to the input frequency and damping 
parameters. The results of this investigation indicate the following conclusions: 

1 . The criterion is simple in concept and application, and it serves as a good 

flutter-onset predictor because it gives a nearly linear variation with 
dynamic pressure. 

2. Emphasis should be given to the accurate measurement of the frequency and 

damping of both modes involved in the flutter condition for reliable 
flutter-onset prediction. 

3. The simplified version using only frequencies predicted a highly noncon- 

servative flutter onset for one case and, thus, should not be used in 
general. 
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4. Normalizing the flutter-margin parameter by the wind-off values of the 

simplified flutter-margin parameter yields a parameter of the order of 
1.0 to 0, and this is the recommended procedure. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 28, 1982 
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Figure 1 


- Geometry of wind-tunnel model. All dimensions are given in inches unless 

otherwise specified. 
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(b) Root locus. 
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(b) Root locus. 



(c) Calculated flutter margin. 


Figure 3.- Concluded 




1.2 

1.0 

.8 

F .6 

.4 

.2 



(a) Variation with . f 1 = 5.23 Hz; = 19.13 Hz. 

Figure 6.- Variation of F with frequency and damping parameters. 
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Figure 10.- Comparison of flutter-margin parameters by using 
calculated and measured first-mode data. M = 0.9. 
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